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Abstract. We study strong-gravity effects on modulation of radiation emerging from accreting compact objects as a possible 
mechanism for flux modulation in QPOs. We construct a toy model of an oscillating torus in the slender approximation 
assuming thermal bremsstrahlung for the intrinsic emissivity of the medium and we compute observed (predicted) radiation 
signal including contribution of indirect (higher-order) images and caustics in the Schwarzschild spacetime. We show that 
the simplest oscillation mode in an accretion flow, axisymmetric up-and-down motion at the meridional epicyclic frequency, 
may be directly observable when it occurs in the inner parts of accretion flow around black holes. Together with the second 
oscillation mode, an in-and-out motion at the radial epicyclic frequency, it may then be responsible for the high-frequency 
modulations of the X-ray flux observed at two distinct frequencies (twin HF-QPOs) in micro-quasars. 
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1. Introduction 

X-ray radiation coming from accreting black hole binary 
sources can show quasi-periodic modulations at two distinct 
high frequencies (> 30 Hz), which appear in the 3 : 2 ratio 
(McClintock & Remillard, 2005). Observations show that the 
solely presence of a thin accretion disk is not sufficient to 
produce these HFQPO modulations, because they are exclu- 
sively connected to the spectral state, where the energy spec- 
trum is dominated by a steep power law with some weak 
thermal disk component. We have shown recently (Bursa et 
al., 2004) that significant temporal variations in the observed 
flux can be accomplished by oscillations in the geometrically 
thick flows, fluid tori, even if they are axially symmetric. Here 
we propose that the QPO variations in the energetic part of 
the spectrum may come from such very hot and optically thin 
torus terminating the accretion flow, which exhibits two basic 
oscillating modes. 

Relativistic tori will generally oscillate in a mixture of in- 
ternal and global modes. Internal modes cause oscillations of 
the pressure and density profiles within the torus. The out- 
going flux is therefore directly modulated by changes in the 
thermodynamical properties of the gas, while the shape of the 
torus is nearly unchanged, which is off our interest here. Glo- 
bal modes, on the other hand, alter mainly the spatial distribu- 
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tion of the material. Because light rays do not follow straight 
lines in a curved spacetime, these changes can be displayed 
out by effects of gravitational lensing and light bending. 

In this paper we summarize extended results of numerical 
calculations and show how simple global oscillation modes of 
a gaseous torus affect the outgoing flux received by an static 
distant observer in the asymptotically flat spacetime and how 
the flux modulation depends on the geometry and various pa- 
rameters of the torus. In Section 2 we briefly summarise the 
idea of the slender torus model and the equations, which are 
used to construct the torus and to set its radiative properties. 
In Section 3 we let the torus to execute global oscillations 
and using a numerical ray-tracing we inspect how these os- 
cillations modulate the observed flux. If not stated otherwise, 
we use geometrical units c = G = 1 throughout this paper. 

2. Slender torus model 

The idea of a slender torus was initially invented by Madej & 
Paczynski (1977) in their model of accretion disk of U Gemi- 
norum. They noticed that in the slender limit (i.e. when the 
torus is small as compared with its distance) and in the New- 
tonian potential, the equipotential surfaces are concentric cir- 
cles. This additional symmetry induced by a Newtonian po- 
tential allowed Blaes (1985) to find a complete set of normal 
mode solutions for the linear perturbations of polytropic tori 



©2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



850 



Astron. Nachr. / AN 326, No. 9 (2005) / www.an-joumal.org 



with constant specific angular momentum. He extended cal- 
culations done for a 'thin isothermal ring' by Papaloizou & 
Pringle (1984) and showed how to find eigenfunctions and 
eigenfrequencies of all internal modes. 

Abramowicz et al. (2005) have recently considered glo- 
bal modes of a slender torus and showed that between pos- 
sible solutions of the relativistic Papaloizou-Pringle equation 
there exist also rigid and axisymmetric (m = 0) modes. These 
modes represent the simplest global and always-present os- 
cillations in an accretion flow, axisymmetric up-down and in- 
out motion at the meridional and radial epicyclic frequencies. 

Metric 

Most, if not all, of stellar and super-massive black holes have 
considerable amount of angular momentum, so that the Kerr 
metric has to be used to accurately describe their exterior 
spacetime. However, here we intend to study the basic ef- 
fects of general relativity on the appearance of a moving ax- 
isymmetric body. We are mainly interested in how the light 
bending and gravitational lensing can modulate observed flux 
from sources. For this purpose we are pressing for a maxi- 
mum simplicity to be able to isolate and recognise the essen- 
tial effects of strong gravity on light. 

Therefore, instead of the appropriate Kerr metric, we 
make use of the static Schwarzschild metric for calculations 
and where we compare with the non-relativistic case, the 
Minkowski flat spacetime metric is also used. 



Equipotential structure 

The equipotential structure of a real torus is given by the Eu- 
ler equation. 



(1) 



where a^^u^ViyU^ is the 4-acceleration of the fluid and e, 
p are respectively the proper energy density and the isotropic 
pressure. The fluid rotates in the azimuthal direction with the 
angular velocity Vl and has the 4-velocity of the form 



= {u\ 0, 0, u*") = -u* (1, 0, 0, n) 

After the substitution of (2), the Euler equation reads 



P 



(2) 



(3) 



where ZY = In (g** + i'^g'^'^) is the effective potential and 
£ is the specific angular momentum. 

For a barotropic fluid, i.e. the fluid described by a one- 
parametric equation of state p=p(e), the surfaces of constant 
pressure and constant total energy density coincide and it is 
possible to find a potential W such that W — ~ J^^ Vp/ [p + 
e), which simplifies the problem enormously (Abramowicz, 
Jaroszyriski & Sikora, 1978). The shape of the 'equipotential' 
surfaces W{r^ z) = const is then given by specification of the 
rotation law I — and of the gravitational field. 
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Fig. 1. An illustration of the equipotential structure of a real rela- 
tivistic torus {lower part) and of our circular slender torus model 
(upper part) surrounding a black hole. The equipotential contours 
are separated by equal steps in the potential W . 



We assume the fluid to have uniform specific angular mo- 
mentum. 



ro - 2M 







(4) 



where tq represents the centre of the torus. At this point, grav- 
itational and centrifugal forces are just balanced and the fluid 
moves freely with the rotational velocity and the specific an- 
gular momentum having their Keplerian values f^Kl^'o) and 

The shape of the torus is given by the solution of equa- 
tion (3), which in the case of constant £ has a simple form. 



W = U + const 



(5) 



In the slender approximation, the solution can be expressed 
in terms of second derivatives of the effective potential and 
it turns out that the torus has an elliptical cross-section with 
semi-axes in the ratio of epicyclic frequencies (Abramowicz 
et al. 2005; see also [Sramkova] in this proceedings). 

In the model used here, we make even greater simplifica- 
tion. Introducing the cylindrical coordinates (i, r, z, (f>), we 
use only the expansion at r = ro in the z-direction to obtain 
a slender torus with a circular cross-section of equipotential 
surfaces. 



W{r, z) = i In 



(ro - 2Mf 



ro (ro - 3Af) 



M[ir- 



2rl {ro~3M) 



(6) 



The profiles of the equipotential structure of a relativistic 
torus and of our model are illustrated in Fig. 1. 

Thermodynamics 

An equation of state of polytropic type, 

p = Kp\ (7) 

is assumed to complete the thermodynamical description of 
the fluid. Here, 7 is the adiabatic index, which have a value of 
^/ 3 for an adiabatic mono-atomic gas, and K is the polytropic 
constant determining the specific adiabatic process. 

Now, we can directly integrate the right-hand side of the 
Euler equation (1) and obtain an expression for the potential 
W in terms of fluid density. 



ly = In p - In (isT 7 p''' + p 7 - p) + In (7 - 1) 



(8) 
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where we have fixed the integration constant by the require- 
ment W{p — 0) — 0. The density and temperature profiles are 
therefore 



P 



T : 



7-1 



-1) 



mu /^w 7 



1 



1 



(9) 



(10) 



ks p A:b 7 
where ks and TOu and the molecular weight, the Boltz- 
mann constant and the atomic mass unit, respectively (Fig. 2). 

Bremsstrahlung cooling ' 

We assume the torus to be filled with an optically thin gas ra- 
diating by the bremsstrahlung cooUng. The emission include 
radiation from both electron-ion and electron-electron colli- 
sions (Stepney & Guilbert, 1983; Narayan & Yi, 1995); 

f = fe.+fee. (11) 

The contributions of either types are given by 
fei = UefiaT caj meC^ Fei{9e) and (12) 
fee = nlcrjaf rUe C^Fee{0e) , (13) 
where rig and n are number densities of electrons and ions, 
is Thomson cross-section, nie and Ve — e^/rrieC^ denote mass 
of electron and its classical radius, is the fine structure 
constant. Fee (fie) and F^i {0^ ) are radiation rate functions and 
Oe — kTe/meC^ is the dimensionless electron temperature. 
Fee{Oe) and Fei{Oe) are about of the same order, so that the 
ratio of electron-ion and electron-electron bremsstrahlung is 

1^ « ^ « 8.4 (14) 

Jee fe 

and we can neglect the contribution from electron-electron 
collisions. For the function F^i (Be) Narayan & Yi (1995) give 
the following expression: 

= 4 f^y^' [1 + 1.78101-34] , 9e<l, 



(15) 



9^ 
2n 



[ln(1.123 6ie + 0.48) + 1.5] , 0e > I 



(16) 



CGS units are used in this paragraph 
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Fig. 2. The density (left) and temperature (right) profiles of a poly- 
tropic gas forming an accretion torus with the centre at ro = 10.8 M. 
Solid lines represent the slender model with radius Ro = 2 M and 
dashed lines represent the real torus filling the potential well of the 
same depth. 




Fig. 4. A schematic illustration of the displacement. The centre T of 
the torus is displaced radially by 5r and vertically by 5z from its 
equilibrium position E, which is at the distance ro from the centre 
of gravity G. 



In case of a multi-component plasma, the density n is cal- 
culated as a sum over individual ion species, n = J2 ^j' 
where Zj is the charge of j-th species and Uj is its num- 
ber density. For a hydrogen^elium composition with abun- 
dances X : Y holds the following: 



Ue = Zj Uj 



Ai-j mu Uj 



X+2Y 
X+Y 


E 


nj 


X+4Y 
X+Y 


E 




X+4Y 
X+Y 


E 





(17) 
(18) 
(19) 



where A^j is the relative atomic weight of the j-th spe- 
cies, TOu denotes the atomic mass unit and we define 
p, = {X + AY) /{X + Y). The emissivity is then 



/e, = 4.30 X 10-25 p2 p^^^Q^^ gj.g 



(20) 



which for the non-relativistic limit {0,. <^ 1) and Population I 
abundances (X — ().l and Y — 0.28) gives 



fei = 3.93 X 10^" p2 T'^ erg cm"'"^ s^^ 



3. Torus oscillations 



(21) 



In the following, we excite in the torus rigid and axisymmet- 
ric (m — 0) sinusoidal oscillations in the vertical direction, i.e. 
parallel to its axis, as well as in the perpendicular radial direc- 
tion. Such an assumption will serve us to model the possible 
basic global modes found by Abramowicz et al. (2005). In 
our model, the torus is rigidly displaced from its equilibrium 
(Fig. 4), so that the position of the central circle varies as 



r(t) — 7'o + 5r sin(wrt) , z(t) — Sz sm{u!zt) 



(22) 



Here, uj^ = = (M/rg)^ is the vertical epicyclic fre- 
quency, in Schwarzschild geometry equal to the Keplerian 
orbital frequency, and ujr ~ ^k{^ — QM/ro)^ is the ra- 
dial epicyclic frequency. The torus is placed at the distance 
ro = 10.8 M so that the oscillation frequency ratio ujz : Wr is 
3 : 2, but the choice is arbitrary. If not stated otherwise, the 
cross-section radius is i?o = 2.0 M and amplitudes of the both 
vertical and radial motion are set to Sz — 5r ^ 0.1 Rq. 

We initially assume the 'incompressible' mode, where the 
equipotential structure and the thermodynamical quantities 
describing the torus are fixed and do not vary in time as the 
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Frequency [Hz] Piequeiicy [Hz] Frequency [Hz| 

Fig. 3. Power spectra of an oscillating torus calculated in the Newtonian limit (left), Minkowski spacetime {middle) and the Schwarzschild 
spacetime (right). Viewing angle is 70°. 



torus moves. Later in this Section we describe also the 'com- 
pressible' mode and discuss how changes in the torus proper- 
ties affect powers in the different oscillations. 

The radial motion results in a periodic change of volume 
of the torus. Because the optically thin torus is assumed to 
be filled with a polytropic gas radiating by bremsstrahlung 
cooling and we fix the density and temperature profiles, there 
is a corresponding change of luminosity Lex j f dV, with a 
clear periodicity at 2t\ jiOr. On the contrary, the vertical mo- 
tion does not change the properties of the torus or its over- 
all luminosity. We find that in spite of this, and although the 
torus is perfectly axisymmetric, the flux observed at infinity 
clearly varies at the oscillation frequency a;^. This is caused 
by relativistic effects at the source (lensing, beaming and time 
delay), and no other cause need to be invoked to explain in 
principle the highest-frequency modulation of X-rays in lu- 
minous black-hole binary sources. 

Effect of spacetime geometry 

In the Newtonian limit and when the speed of light c-^oo, the 
only observable periodicity is the radial oscillation. There is 
no sign of the lo^ frequency in the power spectrum, although 
the torus is moving vertically. This is clear and easy to un- 
derstand, because the c^oo limit suppresses the time delay 
effects and causes photons from all parts of the torus to reach 
an observer at the same instant of time, so it is really seen as 
rigidly moving up and down giving no reason for modulation 
at the vertical frequency. 

When the condition of the infinite light speed is relaxed, 
the torus is no longer seen as a rigid body. The delay between 
photons, which originate at the opposite sides of the torus at 
the same coordinate time, is Ai ~ 2 tq/c sin z, where i is the 
viewing angle {i.e. inclination of the observer). It is maximal 
for an edge-on view V2) ™d compared to the Keplerian 
orbital period it is A</Tk ~ (27r^ ro/rg)^^/^. It makes about 
10% at ro = 10.8A/. The torus is seen from distance as an 
elastic ring, which modulates its brightness also at the vertical 
oscillation frequency Wz due to the time delay effect and the 
seeming volume change. 

Curved spacetime adds the effect of light bending. Pho- 
tons are focused by the central mass's gravity, which leads to 
a magnification of any vertical movement. Black hole is not 
a perfect lens, so that the parallel rays do not cross in a single 
point, but rather form a narrow focal furrow behind it. When 
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Fig. 5. The inclination dependence of powers in the radial (red) 
and the vertical (blue) oscillations. Top panel shows calculations in 
the flat spacetime, bottom panel shows powers as computed in the 
curved Schwarzschild spacetime. Dashed lines represent the same 
calculations done with switched-off g-factor (g = 1). 



the torus trench the fuiTow (at high viewing angles), its os- 
cillations are greatly amplified by the lensing effect. This is 
especially significant in the case of the vertical oscillation, as 
the bright centre of the torus periodically passes through the 
focal line. 

Figure 3 illustrates the geometry effect on three Fourier 
power density spectra of an oscillating torus. The spectra are 
calculated for the same parameters and only the metric is 
changed. The appearance of the vertical oscillation peak in 
the 'finite light speed' case and its power amplification in the 
relativistic case are clearly visible. 



Effect of inclination 

In previous paragraphs we have find out that both the time 
delay and the lensing effects are most pronounced when the 
viewing angle is rather high. Now we will show how much 
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Torus size i?„ [GA//c-] 



Fig. 6. Powers in the radial {top) and vertical {middle) oscillations 
and their ratio {bottom) as a function of the torus size. Different 
viewing angles are plotted. 




Torus distauce Tq [GA1/c^] 

Fig. 7. Powers in the radial {top) and vertical {middle) oscillations 
and their ratio {bottom) as a function of the torus distance from the 
gravity centre. Different viewing angles are plotted. 



is the observed flux modulated when the torus is seen from 
different directions. 

The effect of inclination is probably the most featured, in 
spite of it is difficult to be directly observed. Changing the 
line of sight mixes powers in amplitudes, because different 
effects are important at different angles. When the torus is 
viewed face-on {i.e. from the top), we expect the amplitude 
of LOr to be dominant, as the radial pulsations of the torus 
can be nicely seen and light rays passing through the gas are 
not yet strongly bended. When viewed almost edge-on, the 
Doppler effect dumps the power of Ur and gravitational lens- 
ing amplifies the power in lOz- Thus we expect the vertical 
oscillation to overpower the radial one. 

Figure 5 shows the inclination dependence of oscilla- 
tion powers in the flat Minkowski spacetime {top) and in the 
curved Schwarzschild spacetime {bottom). We see that in the 
flat spacetime the power of radial oscillation gradually de- 
creases, which is caused by the Doppler effect {c.f. the red 
dotted line in the graph). The vertical oscillation decreases 
as well, but it is independent on the g-factor. At inclinations 
i > 75° it has a significant excess caused by the obscuration 
of part of the torus behind an opaque sphere of radius 2M 
representing the central black hole. 

When gravity is added, the situation at low inclinations 
(up to i~25°) is very much similar to the Minkowski case. 
The power of gravitational lensing is clearly visible from the 
blue line, i.e. the vertical oscillation, progression. It is rais- 
ing slowly for inclinations i>45°, then it shows a steeper 
increase for i > 75°, reaches its maximum at i = 85° and it fi- 
nally drops down to zero. At the maximum it overpowers the 



radial oscillation by a factor of 40, while it is 20 x weaker if 
the torus is viewed face-on. The rapid decrease at the end is 
caused by the equatorial plane symmetry. If the line of sight 
is in the = ^/2 plane, the situation is the same above and be- 
low the plane, thus the periodicity is 2 Wz. The power in the 
base frequency drops abruptly and moves to overtones. 

Effect of the torus size 

The effect of the size of the torus is very important to study, 
because it can be directly tested against observational data. 
Other free model parameters tend to be fixed for a given 
source {e.g. like inclination), but the torus size may well vary 
for a single source as a response to temporal changes in the 
accretion rate. 

The power in the radial oscillation is coiTelated with its 
amplitude, which is set to Sr — 0.1 Rq and grows with the 
torus size. It is therefore evident, that the radial power will 
be proportional to i?o squared. If the amplitude was constant 
or at least independent of Rq, the oj^ power would be inde- 
pendent of i?o too. Thus the non-trivial part of the torus size 
dependence will be incurred by vertical movements of the 
torus. 

Figure 6 shows the PSD power profiles of both the ra- 
dial and vertical oscillations for several different inclinations. 
Indeed, the radial power has a quadratic profile and is more 
dominant for lower viewing angles, which follows from the 
previous paragraph. The power in the vertical oscillation is 
at low inclinations also quadratic and similar to the radial 
one, but the reason is different. The time delay effect causes 
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apparent deformations from the circular cross-section as the 
torus moves up and down, i.e. to and from the observer in the 
case of a face-on view. The torus is squeezed along the line 
of sight at the turning points and stretched when passing the 
equatorial plane. Deformations are proportional to its size, 
being the reason for the observed profile. At high inclinations 
the appearance of strong relativistic images boosts the verti- 
cal oscillation power even more. But, as can be clearly seen 
from the 85° line and partially also from the 80° line, there 
is a size threshold, beyond which the oscillation power de- 
creases though the torus still grows. This corresponds to the 
state, where the torus is so big that the relativistic images are 
saturated. Further increase of the torus size only entails an in- 
crease of the total luminosity, while the variability amplitude 
remains about the same, hence leading to the fractional rms 
amplitude downturn. 

Effect of the torus distance 

The distance of the torus also affects the intensity of modula- 
tions in observed lightcurves (Fig. 7). The power in the radial 
oscillation is either increasing or decreasing, depending on 
the inclination. Looking face-on, the g-factor is dominated 
by the redshift component and the power in ujr is increasing 
with the torus distance being less dumped. When the view is 
more inclined, the Doppler component starts to be important 
and the oscillation looses power with the torus distance. The 
critical inclination is about 70°. 

The power of vertical oscillation generally decreases with 
the torus distance. It is made visible mainly by the time de- 
lay effect and because with the increasing distance of the 
torus the oscillation period also increases, the effect is loosing 
on importance. An exception is when the inclination is very 
high. The large portion of visible relativistic images causes 
the vertical power first to increase up to some radius, beyond 
which it then decays. Both small and large tori do not have 
much of visible secondary images, because they are either too 
compact or they are too far The ideal distance is about 11 AI 
- this is the radius, where the torus has the largest portion of 
higher-order images, corresponding to the maximum of the 
vertical power in Fig. 7. 

Generally, the relative power of the vertical oscillation is 
getting weaker as the torus is more and more far-away from 
the graviting centre. This is most significant for higher view- 
ing angles, where the drop between 8M and 16Af can be 
more than one order of magnitude. On the other hand, for low 
inclinations the effect is less dramatic and if viewed face-on 
the power ratio is nearly independent from the distance of the 
fluid ring. 

Effect of radial luminosity variations 

As already mentioned above, the volume of the torus changes 
periodically as the torus moves in and out. In the incompress- 
ible torus, which we have considered so far, this results in 
a corresponding variance of the luminosity, linearly propor- 
tional to the actual distance of the torus r{t) from the centre. 
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Fig. 8. The inclination dependence of powers in the radial (red) and 
the vertical (blue) oscillations in the compressible mode. This is 
the same figure as Fig. 8, except that it is computed with the in- 
clusion of density scaling. Top panel shows calculations in the flat 
spacetime, bottom panel shows powers as computed in the curved 
Schwarzschild spacetime. Dashed lines represent the same calcula- 
tions done with switched-off (/-factor. 



Because we do not change the thermodynamical properties, 
it also means that the total mass M = J p dV contained within 
the torus is not conserved during its radial movements, which 
is the major disadvantage. In this paragraph we relax this con- 
straint and explore the compressible mass conserving mode. 

A compressed torus heats up, which results in an increase 
of its luminosity and size. These two effects go hand-in-hand, 
however to keep things simple we isolate them and only show, 
how powers are affected if we only scale the density and tem- 
perature without changing the torus cross-section. 

We allow the torus to change the pressure and density 
profiles in a way that it will keep its total mass constant. The 
volume element dV is proportional to r, so that in order to 
satisfy this condition, the density must be scaled as 



p(r, z, t) p°(r, z) 



Tit) 



(24) 



where p° refers to the density profile of a steady non- 
oscillating torus with central ring at radius tq . If we substitute 
for the emissivity from (21), we find out that the luminosity 
now goes with r as 



LW^//(p)dy^/pV3dy^r(t) 



-1.33 



(25) 



L{t) - //dV" - r{t) - 5r sin(tj^i) 



(23) 



The negative sign of the exponent causes the luminosity to 
increase when the torus moves in and compresses. Moreover, 
the luminosity variance is stronger than in the incompressible 
case, because of the greater absolute value of the exponent. 

Figure 8 shows the inclination dependence of oscillation 
powers in the compressible case. Compared to Fig. 5 we 
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Fig. 9. The total observed bolometric luminosity of a steady (non- 
oscillating) torus as a function of inclination. In a flat spacetime 
(orange) with only special relativistic effects, the total luminosity 
is increased by a factor of two if the view is changed from face- 
on to edge-on. It is even more in a curved spacetime (blue), where 
the relativistic images make a significant contribution. For compari- 
son also calculations with switched-off g-factor (with g being set to 
unity) are shown (dashed lines). 

see that the signal modulation at vertical frequency is not 
affected, but the slope of the radial oscillation power is re- 
versed. A key role in this reversing plays the (/-factor, which 
combines effects of the Doppler boosting and the gravitation 
redshift. 

The Doppler effect brightens up the part of the torus, 
where the gas moves towards the observer, and darkens the 
receding part. This effect is maximal for incUnations ap- 
proaching i-e. for edge-on view. On average, i.e. in- 
tegrated over the torus volume, the brighten part wins and 
the torus appears more luminous when viewed edge-on (see 
Fig. 9). 

The redshift effect adds the dependence on the radial dis- 
tance from the centre of gravity, which is an important fact 
to explain the qualitative difference between Figs. 5 and 8. 
In the incompressible mode, the luminosity has a minimum 
when the torus moves in and a maximum when it moves out 
of its equilibrium position. The g-factor goes the same way 
and consequently amplifies the amplitude of the luminosity 
variability. The situation is right opposite in the compress- 
ible mode and the luminosity has a maximum when the torus 
moves in and a minimum when it moves out. The g-factor 
goes with the opposite phase and dumps the luminosity am- 
plitude. Because the difference in the g-factor value is more 
pronounced with inclination, it results in increasing or de- 
creasing dependence of the radial power on inclination in the 
compressible or incompressible case, respectively. 

4. Discussion and Conclusions 

We have found out that intrinsic variations of the radiation 
emitted from inner parts of an accretion flow may be sig- 
nificantly modified by effects of a strong gravitational field. 
Above all we have shown that orientation of the system with 
respect to the observer is an important factor, which may al- 
ter the distribution of powers in different modes. However 



this effect, although strong, cannot be directly observed, be- 
cause the inclination of a given source is fixed and mostly 
uncertain. 

Within the model there are other parameters, which may 
be used for predictions of powers in different frequencies. 
We have shown that the size of the torus affects the power of 
the vertical oscillation. In this model this corresponds to an 
emission of harder photons from a hotter torus and provides 
a link between the model and observations. From those we 
know (Remillard et al., 2002) that the higher HFQPO peak 
is usually more powerful than the lower one in harder spec- 
tral states, which is consistent with the model, but the exact 
correlation depends on amplitudes of both oscillations. 

The power in the radial oscillation very much depends on 
the thermodynamical properties of the torus and on its be- 
haviour under the influence of radial movements. We have 
shown that different parametrizations of intrinsic luminosity 
in the in-and-out motion (i.e. compressible and incompress- 
ible modes) change power of the radial oscillation. On the 
other hand, the power of the vertical oscillation remains un- 
affected. This is an important fact and it means that the flux 
modulation at the vertical frequency is independent on the 
torus properties, driven by relativistic effects only. 

Another model parameter is the distance of the thin ac- 
cretion disk. The Shakura-Sunyaev disk is optically thick 
and blocks propagation of photons, which cross the equa- 
torial plane at radii beyond its moving inner edge. Most of 
the stopped photons are strongly lensed and carry informa- 
tion predominantly about the vertical mode, thus the pres- 
ence or not-presence of an opaque disk may be important for 
the power distribution in QPO modes. However, this effect 
is beyond the scope of this article and will be described in a 
separate paper. 
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